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LINE-OF-SIGHT STATISTICAL METHODS FOR TURBULENT MEDIUM: 
VCS FOR EMISSION AND ABSORPTION LINES 
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D. Pogosyan, 1 A. Lazarian 2 
RESUMEN 
ABSTRACT 

We present an overview of the Velocity Coordinate Spectrum (VCS), a new technique for studying astrophysical 
turbulence that utilizes the line-of-sight statistics of Doppler-broadened spectral lines. We consider the retrieval 
of turbulence spectra from emission intensity observations of both high and low spatial resolution and find that 
the VCS allows one to study turbulence even when the emitting turbulent volume is not spatially resolved. 
This opens interesting prospects for using the technique for extragalactic research. VCS developed for spectral 
emission lines is applicable to absorption lines as well if the optical depth is used instead of intensity. VCS for 
absorption lines in point-source spectra benefit from effectively narrow beam and does not require dense sky 
coverage by sampling directions. Even strongly saturated absorption lines still carry the information about 
the small scale turbulence, albeit limited to the wings of a line. Combining different absorption lines one can 
develop tomography of the turbulence in the interstellar gas in all its complexity. 
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1. INTRODUCTION 

Turbulence is a key element of the dynamics of 
astrophysical fluids, including those of interstellar 
medium (ISM) , clusters of galaxies and circumstellar 
regions. The realization of the importance of turbu- 
lence induces sweeping changes, for instance, in the 
paradigm of ISM. It became clear, for instance, that 
turbulence affects substantially star formation, mix- 
ing of gas, transfer of heat. Observationally it is 
known that the ISM is turbulent on scales ra nging 
from AUs to kpc (see lElmegreen fc Scaloll2004l) . with 
an embedded magnetic field that influences almost 
all of its properties. 

Using a statistical description is a nearly in- 
dispensable strategy when dealing with turbulence. 
One of the most widely used statistical measures is 
the turbulence spectrum, which describes the distri- 
bution of turbulent fluctuations over scales. For in- 
stance, the famous Kolmogorov model of incompress- 
ible turbulence predicts that the difference in veloci- 
ties at different points in turbulent fluid increases on 
average with the separation between points as a cu- 
bic root of the separation, i.e. \Sv\ ~ I 1 / 3 . In terms of 
direction-averaged energy spectrum this gives the fa- 
mous Kolmogorov scaling E(k) ~ 47rfc 2 P(k) ~ fc 5 / 3 , 
where P(k) is a 3D energy spectrum defined as the 
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Fourier transform of the correlation function of ve- 
locity fluctuations £„(r) = (Sv(x)Sv(x + r)). Other 
types of turbulence, i.e. the turbulence of non-linear 
waves or the turbulence of shocks, are characterized 
by different power laws and therefore can be distin- 
guished from the Kolmogorov turbulence of incom- 
pressible eddies. Substantial advances in our under- 
standing of the scal ing of compressible M HD turbu- 
lence (see review by lCho fc Lazarianl EoOS. and refer- 
ences therein) allows us to provide a direct compar- 
ison of theoretical expectations with observations. 

Recovering the velocity spectra from observations 
is a challenging problem. Observations provide in- 
tegrals of either emissivities or opacities, both pro- 
portional to the local densities, in Position-Position 
Velocity space (henceforth PPV), i.e as a function of 
the direction on the sky and frequency that can be 
viewed as the position in velocity coordinate along 
the line of sight. Statistical measurements in this 
space are affected both by inhomogeneous spatial 
distribution and the motion of the matter. Turbu- 
lence is associated with fluctuating velocities that 
cause fluctuations in the Doppler shifts of emis- 
sion and absorption lines. If turbulence is super- 
sonic, the recovering of its velocity spectrum is pos- 
sible with the Velocity Chan nel Analysis (VCA) of 
the emissio n intensity maps ( Lazarian fc PogosvanI 
20001 12004 LP00,LP04). This technique has been 
successfully applied in several papers starting with 
Stanimirovic fc Lazarian! ( 2001 ). However, the VCA 
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is just one way to use the general description of fluc- 
tuations in the PPV space. One can also study fluc- 
tuations along the velocity coordinate. The corre- 
sponding technique was term e d the Velocity Coor- 
dinate Spect rum in Lazarianl (12004) and has been 
developed in lLazarian fc Pogosvanl (12 006. LP06) for 
emission and in Lazarian &: Pogosvanl (|2008l LP08) 
for absorption lines. In this presentation we give an 
overview of these developments. 

2. LINKING OBSERVATIONAL AND 
TURBULENT STATISTICAL MEASURES 

Let us consider a volume of turbulent gas or 
plasma (a "cloud") which extent along the line of 
sight S is much smaller than the distance to the ob- 
server. Along a line of sight that is labeled by a two- 
dimensional position vector X on the cloud image, 
the intensity in a spectral (either emission or absorp- 
tion) line with rest-frame frequency v$, Ix.{v,vq), is 
determined by the density of emitting (or absorbing) 
atoms in PPV space (LP06,LP08), p s (X,v) 
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where e is the emissivity and a is the absorption 
coefficient, and the velocity v w -^-(v — vq), if the 
intrinsic line width is neglected. 

In the simplest case of optically thin emission 
lines and unsaturated absorption lines the PPV den- 
sity is directly connected to observables, the emission 
intensity or the optical depth for absorption 

p s (~K,v) oc ix(w), a — > 0, emission (3) 
p s (X, v) oc Tx_(v) = — ln(I/Io), absorption(A) 

PPV density p s (X.,v) encodes the information 
about both the turbulent velocities and the inhomo- 
geneity in the matter distribution. The characteris- 
tic feature that allows us to distinguish between the 
thermal and turbulent velocity is the spatially corre- 
lated nature of the latter. Under set of assumptions 
employed in LP00 and LP04, the PPV structure 
function d s (R,v) = ([p s (Xi, v x ) - p s (X 2 , v 2 )] ), 
R = |Xi — Xi|, v — \v\ — i)2 1, can be expressed 
via real space velocity and density correlations as 

f s €(r) r v 2 

d s (R,v) oc / dz — = exp — , N — — 

(5) 



where the structure function of the velocity z- 
components is 



D z (r) a CV" 



(6) 



( C is a normalization constant ) and for the density 
correlation function £(r = |xi — xa|) = (/?(xi)/7(x2)) 
the isotropic power-law form is adopted 



£{r) oc 1 + 



(7) 



The thermal broadening is described by — 

The Eq. ((5|) provides way for recovery of the 
turbulence by measuring the correlations of the ob- 
served signal, for example the structure functions 



V{R,v) = <[/ x >i 
£>(#,<;) = ([r Xl M-Tx 2 K>n 



2 } oc d s {R,v) 



(8) 

Despite formal similarity, use of emission and ab- 
sorption lines has also several distinct features. 

Measuring the emission by the turbulent gas typ- 
ically provides us with a sky map of the signal. 
If the angular resolution of such map is high, it 
is more advantageous to measure angular correla- 
tions across the image V(R, 0). This leads to the 
VCA analysis (LP00, LP04) that has been reviewed 
at the "Magnetic Fields in the Universe I" meeting 
( Pogosvan fc Lazarian| [2005) . The use of the line-of- 
sight measure T>(0,v) that we review here is more 
challenging but has an advantage for the data with 
poor angular resolution. 

Absorption spectra, on the other hand, typically 
represents real pencil beams for which the line-of- 
sight statistics is the only one available. The effec- 
tive angular resolution of such spectra is high if the 
source is a point object. Therefore both high and 
low angular resolution limits of the VCS arc of in- 
terest. One can also applying VCS to the wings of 
saturated absorption lines (LP08) although the is- 
sues of noise, velocity resolution and intrinsic line 
broadening become critical for such studies. 

3. VCS FOR A TRANSPARENT EMITTING 
MEDIUM 

The intensity in optically thin lines provides di- 
rect information on the density in PPV space accord- 
ing to Eq. (jSJ). Finite angular resolution of the in- 
strument leads to the beam-averaged measurements 



V(0,v) = J dKB (R)d s (R,v) 



(9) 



3 Self-absorption in case of emission line 



where the realistic beam has a finite width, AB. 
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We shell formulate VCS results using the velocity 
coordinate power spectrum 

P(k v ) = J dvV(0,v)exp [~k 2 y/2] (10) 

where k v is the "wave-number" reciprocal to the ve- 
locity v coordinate. With regard to the observation 
scale k v , two other scales define the different regimes 
of the analysis: 

From instrumental perspective it is resolution 
scale is important. To the linear scale AS corre- 
sponds the velocity scale 

V AB = VD Z (S)(AB/S) m , (11) 

equal to the magnitude of turbulent velocities at the 
separation of a size AB. For a measurements at the 
scale k v the beam is effectively narrow when 

kv < Vab' 1 (12) 

while on shorter scales its width is important. Note 
that it is angular resolution that affects the be- 
haviour of line-of-sight statistics. Similarly, for 2D 
VCA analysis in channel maps velocity resolution 
played a role in interpreting 2D statistics of inten- 
sity maps 

From the perspective of the turbulence properties 
the amplitude of the density contribution to VCS 
is encoded in the correlation length rg. This scale 
decides whether density inhomogeneities dominate 
in Eq. Q or the turbulence effects are manifest only 
in random motions of the uniform gas distribution, 
£(r) = 1. The velocity scale that corresponds to r 
is 

V ro = <Jd z (S) (ro/S) m , (13) 

So we can roughly equate the scale of equality of 
the density and velocity effects to ro, i.e., in velocity 
units, 

In Figure [1] we summarize the different scalings 
of VCS. As for the VCA the main difference stems 
from the density being either shallow or steep. If 
the density is shallow i.e. scales as £ ~ r~ 7 , 7 > 

0, which means that the correlations increase with 
the decrease of the scale, then it eventually becomes 
important at sufficiently small velocity differences, 

1. e. at sufficiently large k v . In the opposite case, 
i.e. when 7 < 0, the contributions of density can be 
important only at large velocity separations. 

The left and middle panels of Figure [1] deal with 
the case of shallow density. Velocity is dominant at 
k v < fco, while the density term provides the main 



contribution at k v > fco. The left panel demon- 
strates the case where the scale of transition from 
asymptotic is entirely dominated by velocity to the 
one influenced by spatially resolved velocity and den- 
sity, Vab < V ro . The observed fluctuations arising 
from the unresolved turbulent eddies depends on the 
scalings of both velocity and density. In the middle 
panel, Vab > V ro , the transition scale is unresolved. 
In this case if there is still a dynamical range for mod- 
erately long scales to be resolved by the experiment 
D Z {S)^ 2 > k- 1 > Vab, the VCS of the resolved 
eddies will be determined by the turbulent velocities 
only. 

The right panel of Figure Q] addresses the case 
of a steep density spectrum. The difference now 
is that fluctuations of the density are maximal at 
low wave-numbers and it is there that the density 
could be important. Velocity is dominant at shorter 
scales k v > k . However, the steep density correla- 
tion length r is large, at least as large as the density 
power cutoff r c , which argues for density fluctuations 
to be subdominant everywhere up to the scale of the 
emitting turbulent volume ("cloud"). 

3.1. Thermal effects and inertial range 

In Figure [T] we have plotted theoretical power 
spectra over 3 decades of velocity magnitude to com- 
pactly demonstrate all possible scalings. Thermal ef- 
fects are shown as well, which in a language of the 
power spectrum are simply described by exp[— (3k 2 ] 
cutoff factor. 

In observations, one can anticipate a coverage 
over two decades of velocity magnitudes before ther- 
mal effects get important. Potentially, correcting for 
the thermal broadening, one can extend the obser- 
vational range further. Note, that the thermal cor- 
rections are different for species of different mass. 
Therefore, by using heavier species one can extend 
the high k cut-off by the square root of the ratio of 
the mass of the species to the mass of hydrogen. 

How deeply subsonic scales can still be probed in 
cold gas depends on the signal-to-noise ratio of the 
available data. Indeed, the ability to deconvolve the 
thermal smoothing is limited by noise amplification 
in the process. However, any extension of the VCS 
by a factor a in velocity results in the extension of 
the sampled spatial scales by a factor of a 2 /™, which 
is a 3 for the Kolmogorov turbulence. 

Even with a limited k v coverage, important re- 
sults can be obtained with VCS, especially if observa- 
tions encompass one of the transitional regimes. For 
example, if one measures a transition from a shal- 
low spectrum to a steeper one (see the left and the 
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Fig. 1. Qualitative representation of the density and velocity contributions to the VCS power spectrum and the resulting 
scaling regimes. In every panel light lines show contributions from the p-term (density modified by velocity, dashed line) 
and u-term (pure velocity effect, solid line) separately, while the dark solid line shows the combined total VCS power 
spectrum. Thermal suppression of fluctuations is shown by the dotted line. The labels above the dark solid curve are 
arranged so as to illustrate the sequential transition of the scalings of the total power spectrum. Everywhere except the 
intermediate regimes, the total spectrum is dominated by one of the components to which the current labeled scaling 
corresponds. Labels below the dark solid lines mark the scaling of the subdominant contributions. For the left and 
middle panels the density power spectrum is taken to be shallow, 7 > 0. The left panel corresponds to high amplitude 
of the density correlations, ro > AS, where density effects become dominant at relatively long wavelengths for which 
the beam is narrow. In the middle panel, the amplitude of density correlations is low ro < AB and they dominate only 
the smallest scales which results in the intermediate steepening of the VCS scaling. The right panel corresponds to the 
steep density spectrum. In this case the density contribution is always subdominant. In this example the thermal scale 
is five times shorter than the resolution scale Vab- 



right panels in Figure [T]) one has the potential to i) 
determine the velocity index to, since the difference 
between the slopes is always 4/m; ii) determine 7 
next; iii) estimate the amplitude of turbulent veloc- 
ities from the position of the transition point as dis- 
cussed above. On the other hand, if one encounters a 
transition from a steep to a shallower spectrum, one 
i) may argue for the presence of the shallow density 
inhomogeneities; ii) estimate 7/TO from the differ- 
ence of the slopes, and then to; iii) estimate the den- 
sity correlation radius ro. Finally, if no transition 
regime is available, then for 7 < one can get to, 
while for the case of 7 > a combination of 7 and 
to is available. 

4. VCS FOR ABSORPTION LINES 
Studies of turbulence with absorption lines are 
possible with the VCS technique if, instead of inten- 
sity I{v), one uses the optical depth t(v), i.e., the 
logarithm of the absorbed intensity log/ a 6 S (^). 

In the weak absorption regime, i.e. when the 
optical depth at the middle of the absorption line 
is less than unity, the analysis of the t(v) coincides 
with the analysis of emission discussed above, taken 
in the limit of ideal angular resolution. 

In Figure O we summarize different regimes for 
studying turbulence using absorption lines from 



point sources for which angular resolution can be 
considered as perfect. For shallow density, i.e. for 
£(r) n 2 [1 + (ro/r) 7 ], 7 > the spectrum of 
optical depth fluctuations that arise from inhomo- 
geneities in density projected onto velocity coordi- 
nate scales as P p (k v ) ~ k v 2( * 3 7 ^ m . For sufficiently 
large k v its contribution dominates the contribution 
Pv(k v ) ~ k~ 2 l m coming from purely velocity pro- 
jection effect of turbulcntly moving uniformly dis- 
tributed, £(z) ps ft 2 , gas. If we measure the point 
AV ro re y/DJS)(r /S) m / 2 at which P p becomes 
dominant over P v , we can estimate the density cor- 
relation radius ro which has the physical meaning 
of the scale at which the dispersion of the fluctua- 
tions of density equal the mean density (see LP06). 
However, whether we can observe both regimes or 
just one, depends on the mask that is imposed by 
the saturation of the absorption line, as well as the 
small scale filtering arising from thermal broadening. 
For 7 < 0, P v always dominates. 

In the intermediate absorption regime, i.e. when 
the optical depth at the middle of the absorption line 
is larger than unity, but less than 10 3 , the wings of 
the absorption line can be used for the analysis. The 
central, saturated part of the line is expected to be 
noise dominated. 
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Fig. 2. Power spectrum of the optical depth fluctuations 
due to turbulent motions of absorbers in unsaturated ab- 
sorption line. Thermal and intrinsic line smoothing is 
factored out. Effect of density inhomogeneities becomes 
dominant at k v > V^ 1 , assuming shallow density, 7 > 0. 

Higher the absorption, less the extend of the 
wings available for the analysis is. In terms of the 
mathematical setting this introduces and additional 
window in the expressions for the VCS analysis. 
However, the contrast of the small scale fluctuations 
increases with the decrease of the window. 

For strong absorption regime, the broadening 
is determined by Lorentzian wings of the line and 
therefore no information on turbulence is available. 

5. VCS COOKBOOK 

The VCS cookbook for emission lines is rather 
straightforward. VCS near a scale k v depends on 
whether the instrument resolves the correspondent 
spatial scale [fc^D^S')] S, where S is the scale 
where turbulence saturates. If this scale is resolved 
then P v (k v ) oc ky 2 ^ m and P p (k v ) oc ■ If 

the scale is not resolved then P v (k v ) oc k v 6 and 
P p (k v ) oc k v 2 ' 3 J " m . These results are presented in 
a compact form in Table [1] The transition from the 
low to the high resolution regimes happens as the 
velocity scale under study gets comparable to the 
turbulent velocity at the minimal spatially resolved 
scale. As the change of slope is the velocity-induced 
effect, it is not surprising that the difference in spec- 
tral indexes in the low and high resolution limit is 
4/m for both P v and P p terms, i.e it does not depend 
on the density This allows for separation of the veloc- 
ity and density contributions. For instance, Figure [T] 
illustrates that in the case of shallow density both 
the density and velocity spectra can be obtained. 
Note that obtaining the density spectrum from 



a well resolved map of intensities is trivial for the 
optically thin medium, as the density spectrum is 
directly available from the column densities (i.e. ve- 
locity integrated intensities). However, for the ab- 
sorbing medium such velocity-integrated maps pro- 
vide the universal spectrum K~ 3 , where K is the 
2D wavenumbcr (LP04). Similarly, even for the op- 
tically thin medium, it is not possible to get the den- 
sity spectrum if the turbulent volume is not spatially 
resolved. On the contrary, P p (k v ) reflects the con- 
tribution of shallow density even in this case. 

An advantage characteristic to absorption lines 
is that they provide measurements in the "high 
resolution regime" , provided that the absorption 
against a point source is used. This is impor- 
tant, in particular, for turbulence research in extra- 
galactic setting where even high angular resolution 
may not translate into high spatial resolution. To 
get proper statistical averaging, one requires more 
than one point source. Numerica l experiments in 
Chepurnov &: Lazarian ( 20061 20081) showed that us- 



ing measuring spectra along 5-10 lines of sight is 
enough to recover the underlying spectrum properly. 

In the simplest case, Figure [2] illustrates how the 
information on velocity and density can be obtained 
from observations for the case of shallow density. In 
the case of steep density the measured spectrum is 
always ~ k v 2 ^ m , which allows finding the spectral 
index of velocity m. 

Similar to the case of emitting gas, the VCS study 
of absorbing gas is dominated by the coldest gas 
component. A practical e mission data handl i ng us - 
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ing the VCS technique in iChepurnov et al. 
showed that fitting the observed power spectrum us- 
ing the integral expressions is a more reliable way of 
obtaining both parameters of the underlying turbu- 
lence and temperature of the emitting gas. 

The difference in data analysis for abs o rption 
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lines compared to that of IChepurnov et al. 
is twofold. Firstly, one has to take into account 
masking of the data that restricts it to the wings 
region if the line is sa turated. Secondly , for e mis- 
sion maps analyzed in Chepurnov et al. ( 20061 ) the 
fitting was done while varying the spatial resolution 
of the data. As the effects of particular factors, e.g. 
thermal broadening, are different at different spa- 
tial scales, one can provide additional testing of the 
accuracy of the fit by varying the angular resolu- 
tion. Varying the resolution is not an option for the 
absorption data from point sources but is possible 
when spatially extended sources are used to sample 
the medium. In the latter case the VCA-type studies 
of the absorption lines wings may also be possible. 
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TABLE 1 

SCALINGS OF VCS P{K V ) = P p {k v ) + P v {k v ) FOR SHALLOW AND STEEP DENSITIES. 
Spectral term AB < S [k%D z (S)]~™ AB > S [k 2 v D z (S)]~™ 

/ \-2(l--y)/m / \-2(3-7)/m 

P P (K) oc (k v Dl /2 {S)) oc (k v Dl /2 (S)j 

P v {k v ) ot(k v Di /2 (s)y 2/m ot(k v Di /2 (s)y 6/m 
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